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Abstract 

We investigate the effects of cosmological constant on the characteristic peaks in the iron line 
profile in the black hole accretion disk X-ray spectrum. Our results show that for a fixed mass black 
hole, the peaks become less pronounced and closer together with increasing cosmological constant. 
This effect is mainly due to the slower rotational velocity of Keplerian orbits at large radii in the 
Schwarzschild-de Sitter spacetime as compared to those for the Schwarzschild spacetime. This 
change of the iron line profile is similar to that obtained from increasing the outer radius of the 
accretion disk size or reducing the emission intensity power law exponent in the accretion disk 
emissivity models. 
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I. INTRODUCTION 



Observing the iron K a line in X-ray fluorescence emitted from the accretion disk around 
a black hole is an attractive way to probe the geometry of spacetime in this region, because 
this normally narrow line is seen by a remote observer as a characteristic two-peak profile 
flourescence spectrum . This profile carries the imprint of the spacetime metric in the 
vicinity of the emitting region; the curvature of spacetime affects this line both by Doppler 
broadening from the moving scatterer in the accretion disk and by gravitational redshift of 
the scattered photons. Quantitative calculations of this iron line profile have been previously 
obtained for two natural choices of black hole backgrounds: the Schwarzschild and Kerr 
solution, < 8 ee e.g. Q, g Q. Q). The suiting speet r a a r e datively the sa m e, hut 
are potentially observationally distinguishable with high resolution spectral observations. In 
recent years, with the launch of Chandra and XMM-Newton, the quality of the observed 
spectra is beginning to approach the resolution neccessary to discriminate between these 
two background spacetimes p. 

It is natural to ask whether or not other background spacetimes will produce identifi- 
able features in precision measurements of iron line profiles. In this paper we will derive 
the formalism needed for the calculation of the spectrum of a narrow emission line in a 
static spherically symmetric metric. We then apply this formalism to consider the effects 
of the cosmological constant A on the line shape of the accretion disk spectra. In section 
[H]we present the necessary formalism and in section ITTT1 we apply these results to the case 
of the Schwarzschild-de Sitter metric. In section II VI we descibe our numerical algorithm, 
an adaptation of [?J for numerically computing the spectra. In section |V] we present our 
calculated line profiles. We find that the cosmological constant modifies the line profile for 
values of the dimensionless parameter AM 2 > 10~ 7 that characterizes this effect. At the 
present time, WMAP observations give a value of Q\ = .73 ± .04 8]. This value corresponds 
to A = 1.3 x I0~ 46 km~ 2 in units where G — c — 1. Thus this cosmological constant causes 
deviations from the Schwarzschild line profile only for a black hole that is unrealistically 
massive, 10 19 solar masses. Alternately, for black holes in galactic cores with masses of 10 6 
to 10 8 solar masses, a cosmological constant 23 or more orders of magnitude greater than 
the current WMAP value is needed to cause such deviations. Although these results clearly 
imply that iron line spectra will not provide an independent measurement of the cosmo- 
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logical constant, they do point out that metric deviations from other spherically symmetric 
sources could produce a measurable effect. 



II. THE IRON LINE SPECTRUM 

To calculate the observed X-ray florescence spectrum from an accretion disk, we need to 
add the flux reaching the remote observer from every infinitesimal element of the disk. The 
intensity and the frequency of flux received from each infinitesimal element is determined by a 
variety of factors: the rotational velocity of the emitter, its direction relative to the observer, 
the lensing and the gravitational red shift effects of the black hole, and the emissivity of 
disk at the position of that infinitesimal element. Following e.g. we assume that the 
fluorescing particles travel along the circular timelike geodesies (Keplerian orbits), that the 
disk itself is angularly uniform, stationary and thin, and that the wavelength of the emitted 
photons is short enough that geometric optics applies (See |l| for the discussion of these 
assumptions^) . There are a variety of approaches to explicitly carrying out this calculation: 
Pi, A) Hi) Si- We follow that of Chen and Eardley,0|, with minor modifications as we 
find it intuitive and easily generalized to an arbitrary static spherically symmetric metric. 

The frequency shift g of a photon travelling from emitter to observer is given by 

is -p-u 

9 = — = , (1) 

v e -p ■ u e 

where v e is the frequency of the photon in the emitter's frame, v is the frequency of the 
photon as measured by the observer, p is the 4-momentum of the photon along its null 
geodesic and u Q b s and u em are the 4-velocities of the observer and the emitter respectively. 

The total bolometric flux F seen by an observer is linked to that emitted from the source 
by a factor of f due to conservation of 4-vo,n m e Q. Hence, the observed and emitted 
intensities at a given frequency are related to each other by a factor of g 3 . 

The observed line flux at a given frequency from the accretion disk can be readily ex- 
pressed in terms of the specific intensity r, u e ) in the emitter's frame and the frequency 
shift g: 

F u = [ g 3 I(n,r,u e )dn, (2) 
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FIG. 1: Emission angle a and inclination angle i in the spherical coordinate system of the observer 
(J3J). Here z is the observer's z-axis, z' is the accretion disk's axis, x is the rc-axis in either coordinates 
and i is the disk inclination. 

where the integration is over the angles subtended by the disk image. Note that F v will 
depend on the inclination of the disc relative to the line of sight of the observer. 

To proceed further, explicit expressions for the intensity and frequency shift are necessary; 
we will derive this first for a general static, spherically symmetric metric. These calculations 
can be considerably simplified by a convenient choice of the coordinates. Following [7], we 
use a spherical coordinate system (t, r, 9, 0) centered around the gravitational source and 
z-axis pointing towards the observer. Further, without loss of generality, we choose the 
accretion disk axis to be in the y = plane, with the axis projection towards negative x. 
(See Fig. [TJ based on Fig.l of The inclination angle i is the angle between z-axis and 
the disk axis. 

The general static spherically symmetric metric can be written in these coordinates as 



Note that spherical symmetry implies that all geodesies will be planar. The 4-velocity of a 
circular timelike geodesic in the accretion disk at the radial coordinate r in this metric can 
be expressed using the linear rotational velocity of a circular timelike emitter v 




(3) 




(4) 
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and a factor 7 



7 = vT^P (5) 



as 

1 1 

u a p = 7(t a + -vcosa6 a + -wsinaf), (6) 

r rsmf 

where t a ,r a , # a and a are the coordinate basis vectors and a is the angle between the plane 
of the disk and the plane defined by the emitting point, the center of the disk and the 
observer. It can be expressed in terms of % and as 

cos a = sin i sin 0. (7) 

The polar angle of the intersection in these coordinates is given by 

cot 9 = cos tan z; (8) 

The 4-momentum of a photon emitted from the accretion disk is 



p a = hv(±t° + ^-R( b ~rr« + ^9«), (9) 

where b is the photon impact parameter. The corresponding null geodesic lies in the plane 
connecting the observer, the emitter and the center of the disk. 

There are two natural choices for the 4- velocity u D of a timelike observer in a spherically 
symmetric metric: that corresponding to an observer maintaining a fixed distance d from 
the black hole, 

< = ~T^t a (10) 

VW) 

and that corresponding to a radially freely falling observer at distance d, 



e I , e 2 



where e is a constant determined by the initial radial velocity. Note that if the spherically 
symmetric metric is not asymptotically flat, (fTU|) will not approach the timelike Killing 
vector t a for large d. These two 4-velocities coincide for the choice of zero radial velocity at 
d, e = y/T(d). In the following we will use pijl . 

The frequency shift is conveniently written in terms of two factors: 

9 = 9obs9disk, 

(12) 



where for a given photon impact parameter b, g t, s depends on d and e only, 



9 ° bs = W) - \lw) { W) ~ 1] y 1 _ nd)h ^ (13) 

and the expression for gdisk depends on the emitter's position r only: 

9disk = —r 1 ; — — — -. (14) 

7(1 — -v sin z sin 0) 

We will assume that the iron line emission in the emitter's frame is monochromatic and 
isotropic with no local line broadening. For this case the dependence on radius r of the 
specific intensity is commonly modelled by a power law dependence. Assuming an optically 
thin disk, 

I(n,r,v e ) = -£- l r- p 6(v e -v ). (15) 

The power p is a real number usually taken to be in the range ofl<p<3;/«isa constant 
with units of power. The factor p! is the cosine of the angle between the emitted photon 
and the disk axis in the emitter's frame, 

, gbcosi 

p = r-j- 16 

r sin 9 



An explicit derivation of the form p! in the observer's coordinates is given in Ref. 

The infinitesimal solid angle subtended by the disk is easily expressed in terms of the 
impact parameter b and <p: 

bdbdcj) 

= ( 17 ) 
Combining terms and reexpressing v e — v/g in (|TH|) we have 

F v = -^f g 3r -^r-v5(v - gv Q )bdbdcj> (18) 

47TGT J 0COS2 

using the identity 8{y/g — vq) = gbiv — gu ). The range of integration is over and b 
such that the null geodesies intersect the disk. The coordinates r, 9 of the intersection point 
are implicitly a function of b, <f), and i; a null geodesic of given b and <fi will intersect the 
accretion disc at angle % at a given r and 9. For a given frequency y, contributions to F v 



are from points on the accretion disk with equal value of g (See ll| Fig. 7 for an example 
of such contours for Schwarzschild.). Observe that the intersection of the null geodesic with 
the accretion disc can occur either before or after the turning point. 
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III. THE SCHWARZSCHILD - DE SITTER METRIC 



We now apply the results of section |TT] to the Schwarzschild - de Sitter metric. This 
metric describes a static spherically symmetric vacuum solution of Einstein equations with 
cosmological constant. Cosmologically, this is a good model of an isolated black hole in a 
dark energy-dominated spatially flat universe. The Schwarzschild - de Sitter metric in static 
coordinates is 

ds 2 = (1 - ^ ~ \^)dt 2 + _ J- 2 _ + rW . (19) 

r 3 

The effective potential for this spherically-symmetric metric is 

I 2 M Ml 2 1 4 9 
V "' = 2^ - - " — - 6 Ar (20) 

where I is the conserved orbital angular momentum. Circular orbits correspond to the 

extrema of this effective potential 

M Ar 2 I 2 3M1 2 . . 

^- — --s + — = °- ( 21 ) 

The circular orbit is stable when the extremum is a minimum: 

A., 15M, M ,. 6M, 

4 )<— (1 ). (22) 

For AM 2 = 0, we recover the familiar circular orbit stability limit for the Schwarzschild 
metric: r > 6M [10]. We now examine in more detail the orbit stability for the case 
where the spacetime is close to that of Schwarzschild, namely that of small cosmological 
constant relative to the black hole mass, 2M << a/A/3 and distances far from the black 
hole, d » 2M [l4|. See, for example, [12j for a detailed analysis of geodesies in the 
Schwarzschild-de Sitter spacetime. 

Unlike in the case of the Schwarzschild solution, circular orbits become unstable at 
both small and large radii in the metric. The instability at small radii occurs as for the 
Schwarzshild solution. For small enough AM 2 , the minimum stable circular orbit is ap- 
proximately r min ~ 6M + 648AM 3 . In contrast to the Schwarzshild solution, there is a 
new instability at large radii; the "negative pressure" produced by the cosmological con- 
stant destabilizes the outer circular orbits. From (|22j) we can calculate that the maximum 
possible stable orbit is approximately at 

,3M., i 



"max — ' ' " 



4A 



(23) 



Objects beyond this limit would either find a smaller stable circular orbit or become un- 
bound. This effect limits the value of A for which the Schwarzchild-de Sitter metric makes 
sense as a background for a black hole with an accretion disk. For black holes of between 1 
and 10 10 solar masses, a range from stellar remnant black holes to the mass of active galactic 
nuclei, the current estimated value of A yields an r max of the order of 20 to 40, 000 parsec. 
The linear rotational velocity seen by the observer of an element of the disk at radius r 

is 

[M , 

" = VT-r- (24) 

or slower than the velocity at the same radius in the Schwarzschild metric, v = \J~^r- This 
is as expected, as the cosmological term acts as a repulsive force and partially screens the 
central mass [15j . 

To calculate the spectrum, we need to choose an observer. We will assume that the ob- 
server has zero radial velocity relative to the black hole; with this assumption, the observer's 
4-velocity coincides with (fTUj) and g f, s becomes 

9obs = ; 1 = (25) 
1 _ 2M _ A d2 



d 3 

independent of the integration parameters of (|18|) . The contribution from the disk is 



i 3M 



Qdisk = ; — } ■ (26) 

1 rAr 2 sina 
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As g bs is independent of the integration parameters in (fTSj) it can be factored out of the 
integral; (fTHjl becomes 

F » = 1% / 9Lk^:r-*6(v - gvo)bdbd<p . (27) 



The factor depends on both the distance to the accretion disk d and A; the integrand 
is now independent of d. 



IV. THE NUMERICAL SPECTRUM CALCULATION PROCEDURE 



Each emission point of the accretion disk "seen" by a remote observer corresponds to a 
distinct pair of impact parameter and polar angle (b,<p). Thus sampling the (&,</>) space is 
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FIG. 2: For numerical calculations the accretion disk image is divided into N 2 nearly square 
elements whose boundaries are lines of constant b and constant phi. The area of each element is 
chosen to be proportional to the square of b. 

equivalent to sampling the whole disk. This approach, due to Q, circumvents the need to 
explicitly evaluate the Jacobian of coordinate transformation from the disk coordinates to 
observer coords, In other approaches (Q, g Q) this aacohian is elated either 
numerically or analytically. 

The spectrum calculation algorithm is as follows: 

1. Specify the metric and the disk parameters: the desired number of disk image elements 
to use, the inner and outer disk radii in dimensionless units, ri n ,r out , the inclination 
angle i of the disk and the emission intensity power law exponent p (see equation (jl5|)) . 

2. Sample the (b, 0) space. To accomplish this, we discretize dip = 27r/N and take 
db = 27r b/N. The impact parameter b ranges from zero to a value determined by the 
outer disk radius, chosen so that no null geodesic with greater b will intersect the disk. 
The flux contribution from each element and its frequency shift is taken to be that of 
center of each element. 



9 



3. For each sample point, integrate the null geodesic inward from the observer to to the 
emission point, defined as the point where the geodesic intersects the disk. (Points 
in (6, 0) space with null geodesies found not to intersect the disk at this step are 
discarded.) 

4. Calculate the frequency shift using equation and the spectral flux from the disk 
element using equation (fT5J) . Add the resulting intensity to the appropriate frequency 
point of the spectrum. 

The usual way of numerically integrating the null geodesic is to use the second-order 
geodesic equations in order to avoid complications near turning points, especially for the 
high-deflection angle case (see, for example, Q|, [ll]], [ill). In our case, however, nearly all 
of the flux comes along null geodesies with small to moderate deflection, because the inner 
radius of the accretion disk is at least at 6M. Thus, the null geodesic equation 



can be integrated by first-order methods. To produce a "bounce" from a turning point, we 
change the sign of % from — to + when the expression under the square root in 1251 becomes 
negative. This approach is stable under changes in stepsize and reproduces certain analytic 
deflection results to high accuracy. 

The code implementing this algorithm is available from the authors. 

Note that we modeled the disk image as a collection of small nearly-square elements 
in the space (6,0) (see Fig. |2J), and calculated the flux from each one, assuming uniform 
flux from each piece. The error introduced by this procedure adds noise to the resulting 
line profile. Higher noise level for the same discretization is seen for large outer radii and 
for high inclination angle. This is because, for a given number of disk elements, fewer 
image elements overlap the important inner disk area for a large disk. Similarly, for a large 
inclination angle, fewer elements overlap the disk, causing a higher noise level. However, 
we observed no measurable systematic error even for significant discretization noise. Our 
solution was to increase the number of points (and line profile computation time) from the 
usual choice of 10 6 to 10 8 in such cases. 

For algorithm verification purposes, we calculated line profiles for the Schwarzschild met- 
ric with the parameter sets taken from [| Fig.l, j?J Fig.4a-7a Our simulation results 




(28) 
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for the Schwarzschild case are indistinguishable from those in the original papers. For ex- 
ample, a set of line profile spectra calculated using the parameters taken from Fig. lb of 
is shown on Fig. [3] and replicates the spectra from that paper. We feel confident that our 
implementation of the accretion disk fluorescence spectrum calculations is solid enough to 
apply it to the Schwarzschild - de Sitter case. 



V. RESULTS FOR THE SCHWARZSCHILD - DE SITTER METRIC 

In the Schwarzschild case, the line profile of the black hole normalized to the height and 
position of the blue peak, does not depend on the mass M of the black hole. Instead, it 
depends on the three dimensionless parameters: the ratios of inner and outer accretion disk 
radius to the black hole mass and the emissivity parameter p. The emissivity parameter p 
determines the sensitivity of the line profile to the other disk radii. The line profile for p = 2 
( Q] ) is more sensitive to the choice of inner and outer radii than that for p — 3 ( [3] ) ; this is 
easily understood as the total flux from a fixed radius falls off as r 2 ~ p . 

The line profile in the Schwarzschild-de Sitter metric is similar to that of the Schwarzschild 
case; however the normalized line profile now depends on four dimensionless parameters. The 
additional parameter is AM 2 as discussed in section IIVI In addition, note that the absolute 
position of the blue peak in frequency will in principle depend on both the distance to the 
observer and the cosmological constant. However, for the range of AM 2 in our simulations, 
this effect is negligible. 

Fig. |3] shows the effects of cosmological constant on the normalized accretion disk fluores- 
cence spectra for the disk of inner radius 7M and outer radius 70M, fixed disk inclination 
i = 30° and a constant emissivity index p = 2. The parameter AM 2 is varied from 1 x 10~ 8 
to 2.2 x 10~ 6 . This maximum value of AM 2 corresponds to an maximum outermost stable 
orbit of radius 70M (cf. 12*31) '). This data set is similar to that in Fig.l of 111]. However, for 

no 

most of the data sets used in [1|, lH , the outer radius exceeds the radius of the outermost 
stable circular orbit values of AM 2 > 10~ 7 , which is barely large enough to produce a no- 
ticeable difference in the spectrum. We reduced outer disk radius to just 70M to allow for 
larger AM 2 and more pronounced effects from the cosmological constant. The simulation 
shows that the effects are negligible for AM 2 below about 10 -7 . Above that, the relative 
height of the red peaks get progressively larger with the increase of the cosmological term. 
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In addition, the peaks tend to get closer together. 

Figures 03|7| show the dependence of spectra on the cosmological term for the specified 
values of inclination angles, emissivity index and disk dimensions. The results are similar 
to that illustrated in equation (J24|) for all cases. 

As discussed earlier, the data above are limited to the parameter sets that ensure stable 
circular orbits for all parts of the disk. Much more dramatic differences from Schwarzschild 
spectra can be obtained if the condition of Keplerian orbit stability at the outer disk is 
neglected. Consider, for example, a line profile for AM 2 = 10~ 5 and p = 2 illustrated in 
Fig. |H1 This line profile deviates significantly from the familiar two peak line profile seen 
generally. Of course, it is very likely that extending the accretion disk into region with 
unstable orbits is unphysical. 

VI. DISCUSSION 

One can intuitively understand the effects of the cosmological constant as follows: A 
acts as a negative pressure, reducing the rotational velocity of the disk at a given radius 
in comparison to the Schwarzschild case. The farther away a circular geodesic is from the 
black hole, the more pronounced is the velocity reduction. Now most of the material is 
concentrated in the outer layers of the disk and as the red (blue) peak is produced by the 
parts of the disk moving fastest away (towards) the observer, slower rotational velocity 
farther away means less pronounced peaks. This effect is somewhat similar to that of 
reducing the maximum outer radius or increasing the emissivity index, both of which shift 
the relative contribution to the spectrum in favor of the inner parts of the disk. 

One can estimate the magnitude of the cosmological constant required to see its effect 
on the iron line profile for black holes. Assuming the cosmological constant value consistent 
with the current observations, This value corresponds to A = 1.3 x 10 _46 fcm -2 in units 
where G — c — 1, calculated from the WMAP-based values for the dark energy density 
Q A = 0.73 ± 0.04 and the Hubble constant h = (0.71 ± 0.0A)100km/s/Mpc Q. We find 
that a solar mass black hole, i.e. order of magnitude size of astrophysical black holes with 
observed iron line profiles yields AM 2 = 2.9 x 10 -46 , an undetectable effect. For black 
holes of 10 8 solar masses, roughly the upper limit of galactic black hole masses, yields 
AM 2 = 2.9 x 10~ 30 . To produce a deviation from the Schwarzschild accretion disk iron 
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line profile, this factor has to larger than 10 -7 , which requires either an increase in A by 23 
orders of magnitude or an increase in the black hole mass by over 11 orders of magnitude 
to 3 x 10 19 solar masses. Thus observations of the line shape of the iron line profile will not 
provide sensitive independent bounds on the cosmological constant. This conclusion also 
holds for other homogeneous dark energy models, which, as they have similar dark energy 
densities have similar metrics. 

Although these results clearly imply that iron line spectra will not provide an independent 
measurement of the cosmological constant, they do provide an example how spherically 
symmetric metric deviations can potentially produce a measurable effect. 
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FIG. 3: Spectra calculated for the Schwarzchild metric with parameters taken from Q Fig. lb: 
i = 30°, p = 2, r min = 20M, r max = lOOAf, 200M, 500M, 3000M (left to right). 

[14] There is a limit on how big A can be for an accretion disk to appear at all. For example, for 
large A the term (1 — — |Ar 2 ) is negative for all r, thus making the t coordinate spacelike. 

[15] At r = (^)s the central mass is completely screened, so that a particle can "hover" at that 
distance without rotating at all. This distance, however, is farther from the center than the 
maximum stable circular orbit, and so the equilibrium is unstable. 

Q 

[16] The captions of p| Fig.4 and Fig. 6 appear to be swapped in print. 
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FIG. 4: Spectra calculated for a range of KM for inclination angle i — 30°, p — 2, r m i n — 7M, 
fmax = 70M. Note that the line for AM 2 = 1 x 1CT 7 nearly coincides with that for the Schwarzschild 
case on this and other graphs. 
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FIG. 5: Spectra calculated for a range of AM 2 for inclination angle i — 70°, p — 2, r m i n — 7M, 



r max = 70M. 
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FIG. 6: Spectra calculated for a range of AM 2 for inclination angle i — 30°, p — 1, r m j„ — 7M, 
r max = 70M. 
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FIG. 7: Spectra calculated for a range of AM 2 for inclination angle i — 30°, p — 2, r m j„ — 20M, 
r m „.x = 70M. 



16 



1.2 



AM 2 




FIG. 8: Spectra calculated for AM 2 = 1CT 5 for inclination angle i = 30°, p = 2, r min = 7M, 
r m ax = 70M. Note that the maximum radius exceeds the value allowed by the orbit stability- 
condition for AM 2 = 1 x 1CT 5 . 
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